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Abstract. The notion of quasi-angular momentum is introduced to label the 
eigenstates of a Hamiltonian with a discrete rotational symmetry. This concept is 

recast in an operatorial form where the creation and annihilation operators of a 
Hubbard Hamiltonian carry units of quasi-angular momentum. Using this formalism, 
the ground states of ultracold gases of non-interacting fermions in rotating optical 
lattices are studied as a function of rotation, and transitions between states of different 
quasi-angular momentum are identified. In addition, previous results for strongly- 
interacting bosons are re-examined and compared to the results for non-interacting 
fermions. Quasi-angular momentum can be used to distinguish between these two 
cases. Finally, an experimentally accessible signature of quasi-angular momentum is 
identified in the momentum distributions of single-particle eigenstates. 
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1. Introduction 

Recent experimental and theoretical studies of quantum gases of ultracold atoms have 
been successful in replicating the behavior of a wide class of condensed matter systems. 
In particular, Bose-Einstein condensation (BEC) [H El E] and the degenerate Fermi 
gas |1] have both been observed, and there have been extensive studies of superfluidity, 
vortex formation [5], and the BCS-BEC crossover [6l[7j. The Mott-insulator to superfiuid 
phase transition in the Bose-Hubbard model has been observed in a system of ultracold 
bosons in an optical lattice of [H |9]. In addition, experiments observing novel phases 
of ultracold matter in atomic mixtures of Bose and Fermi gases in optical lattices have 
been performed [lOl [IH HH |13] . 

Along these lines, focus has turned to strongly-correlated effects in systems found 
outside the realm of condensed matter. One promising line of research in the realization 
of these systems is that of coupled arrays of optical cavities, which allow single-site 
addressing [T^ [TB] . Another is a mapping between the Hamiltonian of electrons 
confined in two dimensions in the presence of a constant transverse magnetic field and 
that of a rotating atomic gas. This has led to various theoretical studies predicting 
fractional quantum Hall effect (FQHE) behavior in rotating BEC's [T71 [T8l IT9l 121)]. 
However, experimentally reaching the parameter regimes necessary to observe this 
behavior is difficult [21]. Optical lattices offer a solution to this problem by enhancing 
correlations. Optical lattices are formed as standing waves of counter-propagating laser 
beams and act as a lattice for ultracold atoms. These systems are highly tunable: lattice 
spacing and depth can be varied by tuning the frequency and intensity of the lasers, and 
interactions between atoms can be tuned via a Feshbach resonance |22|. Some recent 
studies have made direct connections between FQHE physics and strongly-interacting 
bosons in optical lattices in the presence of an effective magnetic field [231 1211 [251 ISHl [27] . 

Connecting these theoretical studies to experiment requires identifying observables 
that can act as experimental signatures for the predicted physics. For instance, linear 
response theory has been applied to current fiow in the presence of a "potential gradient" 
in order to observe FQHE physics |26], noise correlation analysis has been applied to 
bosons in a rotating ring lattice [2H], and Bragg spectroscopy has been mentioned as 
a probe for vortex states j29|. A common observable in ultracold atomic gases is the 
momentum distribution, measured in time-of-fiight experiments. 

In this paper, we present analytical results for strongly-interacting bosons in 
rotating optical lattices that confirm numerical results in a previous paper [30| and 
extend them to systems of non-interacting fermions. These results are based on the 
notion of quasi-angular momentum, which is a quantum number for systems with 
a discrete rotational symmetry. Quasi-angular momentum is analogous to quasi- 
momentum for periodic translationally invariant systems |3T] and has been previously 
used in the context of rotating ring lattices [321 [331 131] and carbon nanotubes [35] to 
label eigenstates. A formalism for this quantum number for second-quantized systems is 
presented here and applied to gases of bosons and fermions in rotating optical lattices. 
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The results presented identify transitions between states of different symmetry for the 
ground state of a rotating system. A possible avenue for experimentally detecting these 
results via the momentum distribution of the ground state is also presented. 

The quasi-angular momentum, m, of the ground state of an ultracold quantum gas 
of strongly-interacting bosons or non-interacting fermions in a rotating ring lattice of 
sites can be monitored as a function of rotation speed. For the case of strongly- 
interacting bosons, m cycles through the values m = nl mod : Z = 0, 1, 2, . . ., where 
n is the number of particles. For the case of fermions, m cycles through the values 
m = nl mod A^ : / = 0, 1, 2, . . . for odd numbers of particles and m = n{l + 1/2) mod A^ : 
/ = 0, 1, 2, . . . for even numbers of particles. A system of non-interacting fermions is 
thereby distinguishable from a system of hard-core bosons. Similar behavior obains 
for systems of fermions in two-dimensional square lattices. Signatures of this quantum 
number in experiment are observable via the momentum distribution of the ground 
state. For single-particle systems in two-dimensional square lattice geometries, such 
signatures include the existence of a peak at zero momentum only for the m = state 
and peak-spacing differences between m = 2 and m = 1, 3 state. 

The paper is structured as follows. In section 2, we draw the analogy between 
systems with a discrete translational symmetry and systems with a discrete rotational 
symmetry, thereby generating a Bloch theory for the latter and introducing the notion 
of quasi-angular momentum. Specific investigations of quantum gases of strongly- 
interacting bosons and non-interacting fermions in ring and square lattice geometries 
are carried out in section 3. In section 4, signatures of the quasi-angular momentum 
in the momentum distribution of a state are identified. Section 5 summarizes the main 
results of the paper. 

2. Quasi-angular momentum 

In this section, we briefly review band theory and introduce the language of Bloch 
functions and quasi-momentum (see for instance [361 [37]). By analogy, the quasi-angular 
momentum of an eigenstate of a system with a discrete rotational symmetry is defined. 

2.1. Bloch theory 

When a small periodic potential of lattice period d is introduced to a free-particle 
system, the two wavefunctions ipi {x) = e^'^^l'^ and 1^)2 {x) = e~™^^'^ at the Brillouin zone 
boundaries (g = ivr / d) mix with each other. The potential acts as a perturbation that 
breaks the degeneracy between these two states, creating new eigenfunctions. 



This process opens up a band-gap at the Brillouin zone boundaries. In a reduced zone 
scheme, this energy spectrum becomes the familiar one-dimensional Bloch band diagram 
(see figure 1). 




(1) 
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Figure 1. Dispersion relation for free particles (solid) and particles moving in 
a one-dimensional sinusoidal lattice (dashed) plotted in a reduced zone scheme. 
Free space dispersion relations are plotted versus momentum, and lattice dispersion 
relations are plotted versus quasi-momentum. Energy is in units of the recoil energy, 
E]i — h'^n'^ /2m(P , where d is the lattice spacing. The lattice depth is Vq = 5Ej^. 



The introduction of a periodic potential breaks the continuous translational 
symmetry of the system, and therefore the wavefunctions are no longer momentum 
eigenstates. Instead, the eigenfunctions of the Hamiltonian are eigenfunctions of a 
discrete translation operator, T^, that translates the entire system by one lattice spacing, 

(x - rf) = Td^P (x) = e-"""^ (x) , (2) 

and can be labeled by fc, the quasi-momentum. A quasi-momentum eigenstate of a 
periodic Hamiltonian can be written as the product of a plane wave and a function 
periodic in the lattice: 

^« (x) = e^'^uf (x) , {x-d)= uf (x) , (3) 

where / is a band index. 

When these Bloch functions are expanded in the momentum basis, the only 
momenta contributing to the sum are those that differ from /c by a reciprocal lattice 
vector, Gj = The expansion is then given by 

/oo °° „i(k+2nj/d)x 



where the expansion coefficients c^''' are given by 

"OO -i{k+2nj/d)x 



cf = j Jx "" ^ ^f{x). (5) 

Figure 2(a) is an example of a Bloch function in a sinusoidal lattice, and figure 2(b) is 
the Fourier transform of this Bloch function - that is, the Bloch function in momentum 
space. 

In the rest of the paper, we will consider only lattices that are deep enough that a 
single-band approximation will be adequate to describe the physics. In the proceeding 
sections, then, we drop the band index, /. 
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Figure 2. (Color online.) The zero quasi-momentum Bloch function in the first band 
for an infinite sinusoidal lattice if lattice depth Vq = Si?^. (a) The real space wave 
function is concentrated in the lattice sites, (b) The wave function in momentum space 
exhibits peaks at reciprocal lattice vectors. The envelope of these peaks is given by a 
Gaussian for large lattice depths, consistent with a lattice of simple harmonic oscillator 
wells. 



2.2. Lowering the symmetry 

If an extra potential, V (x), satisfying 

V{x-pd) = V (x) , (6) 

with p > 2 an integer and d the lattice spacing, is introduced to a lattice, then Bloch 
states {x\k) = ipk (x) and {x\k') = ipk' (x) will mix with each other if k and k' differ by 
27rm/pd, as shown by the following calculation: 

/oo 
dx^l, (x) V (x) i)k {x) = e'^^'-^'^P'^ {k' \V\k). (7) 
-oo 

Eigenstates of the new Hamiltonian will therefore be linear combinations of these Bloch 
states. The potential increases the lattice periodicity and, accordingly, decreases the 
size of the first Brillouin zone to 2-^ jpd. 

In figure 3(a), the lowest band of a one-dimensional lattice with an arbitrary 
periodic potential is plotted in a reduced zone scheme. Introducing an infinitesimal 
potential whose periodicity is three times the lattice spacing c?, the first Brillouin zone is 
reduced to the region — 7r/3(i < k <'k /?>d. Equation ([T]) shows that the states lying on a 
vertical line in figure 3(b) mix with each other under the addition of the new potential 
- i.e., {k' \V\ k) is non-zero for these states. In particular, the degenerate states at = 
mix, breaking the degeneracy and opening up a gap at A; = 0, as shown in figure 3(b). 
The single band splits into three bands due to the reduction of the symmetry by V (x). 

When the addition of a potential reduces the symmetry of a system, as V (x) did 
to the generic system described in figure 3, the original eigenstates of the system mix 
according to their symmetry. Accordingly, gaps open up where there are degeneracies 
between states of the same symmetry, creating a band structure. 

2. 3. The ground state of a moving lattice 

The utility of these symmetry concepts and quasi-momentum can be demonstrated 
by considering the effect of moving the lattice. In a stationary lattice, a zero quasi- 
momentum Bloch function is always the ground state. When the lattice is moving. 
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Figure 3. (Color online.) (a) Lowest band dispersion relation for a one-dimensional 
sinusoidal lattice plotted in a reduced zone scheme that anticipates the lowering of the 
symmetry via an additional potential with periodicity 3d. (b) After the introduction 
of this potential, states with the same quasi-momentum, k in the reduced zone mix 
with each other, opening up gaps anywhere that there is a degeneracy. The matrix 
elements of the added potential between states of the same quasi-momentum are of 
magnitude 5i?fl. 



the energetically favored state is one that moves along with the lattice. In order to 
investigate this phenomenon, we move to a frame co-moving with the lattice. The 
Hamiltonian in this frame is given by |38] 

H = Ho-p-^r, (8) 

where Hq is the Hamiltonian in the non-moving frame, and v is the velocity of the 
moving frame with respect to the lab frame. The ground state oi Hq — p ■ \ will be the 
ground state of the lab frame but written in rotating frame coordinates. 

Consider an eigenstate of Hq, ipk (x). In first-order perturbation theory, the energy 
of this state is 

Eu = {k\Ho-p-^^\k) = El-v {k \p\ k) . (9) 

At zero velocity, the ground state is a zero quasi-momentum state, 'ipoi^x). However, 
the energy of a state that has a positive average momentum, {k \p\ k) > 0, falls below 
that of ipo [x) for large enough v. This change is signalled by an exact energy level 
crossing in the energy spectrum between two quasi-momentum states (see figure 4). As 
the velocity is increased, the effect of the velocity term is to reorder the energies of the 
quasi-momentum states, thereby leading to exact energy level crossings. 



2.4- Quasi-angular momentum 

A moving A^-site linear lattice with periodic boundary conditions is analogous to a 
rotating A^-site ring lattice. We can make this analogy explicit by considering the 
Hamiltonian of a moving, one-dimensional, sinusoidal A^-site lattice in the co-moving 
frame, given by 

h"^ d'^ h d 

where v is the velocity of the lattice. The wave- number q can be rewritten as q = n/d, 
where d is the lattice spacing. If periodic boundary conditions, 

ip{x + Nd) = i; (x) , (11) 
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Figure 4. (Color online.) The four lowest energies in the single-particle energy 
spectrum as a function of lattice velocity for a four-site sinusoidal lattice with periodic 
boundary conditions. At zero velocity, the zero quasi-momentum state (solid) is 
the ground state. At some finite positive velocity, the k — Tr/2d state (dashed) is 
energetically favorable and becomes the ground state. The ground state eventually 
cycles through all quasi-momenta. 



are included, then the Hamiltonian describes the system illustrated in figure 5(a) for 
the explicit case of an 8-site lattice. On the other hand, the Hamiltonian of a rotating, 
sinusoidal A^-site ring lattice (figure 5(b)) in the rotating frame is given by 

where is the rotation frequency, R is the radius of the ring, and f ^ is the angular 
mometum operator, L^. The inclusion of the term —QLz has the effect of moving to a 

are 



12 



frame co-rotating with the lattice [38] . The two Hamiltonians, equations 10 and 
identical if we perform a transformation x = (pN (1/271 and identify Nd/27i with R and 
v/R with Q. 

Since the Hamiltonians are exactly identical, all of the properties of one-dimensional 
systems with a discrete translational invariance carry over for ring systems with a 
discrete rotational invariance. The analogy can be carried further, starting with the 
two-dimensional free-space solution in polar coordinates, 

^,(0,p)=e^^*/2,(p), (13) 

where Rj (p) is a radial function, irrelevant for our discussion, and j is an integer. In 
the presence of a potential that breaks the rotational symmetry, the eigenstates are 
linear combinations of these free-space solutions. If this potential has a discrete A^-fold 
rotational symmetry, its eigenstates can be expanded in the free space solutions as 

oo 

i^m{<P,p)= Yl a.e^^^^+-)'^i?, (p) . (14) 

j=-oo 

It is evident from this expansion that ipm {4>, p) is an eigenvector with eigenvalue 
^-t2nm/N q£ ^^^q discretc rotation operator R2-k/n that rotates the system by the 
angle 2tt/N; R2-k/n takes the place of the discrete translation operator, T^. The 
analogy is complete when we note that the eigenstates are linear combinations of 
angular momentum eigenstates, in which case we call the number hm the quasi-angular 
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Figure 5. (Color online.) This diagram shows the equivalence between a one- 
dimensional lattice with periodic boundary conditions and a ring lattice with the same 
number of sites. The translation operator Td that translates through one lattice site 
maps onto the rotation operator i?27r/8 thsA rotates the lattice through one lattice site. 



momentum of the state, ipm {'PiP)- For the rest of the paper, we will drop the factor of 
h from the quasi-angular momentum. 

3. Bosons and fermions in rotating optical lattices 

Interacting bosons in one dimension enter the Tonks-Girardeau regime when interactions 
get very strong [321 SO]- The gas takes on some of the characteristics of a gas of non- 
interacting fermions, such as an identical number density. It differs in certain aspects, 
however, such as the momentum distribution. It then becomes particularly interesting 
to compare the cases of strongly-interacting bosons and non-interacting fermions in 
rotating lattices. In this section, gases of non-interacting fermions and hard-core bosons 
in rotating optical lattices are investigated. A many-particle formalism is generated for 
quasi-angular momentum states and used to monitor the symmetry of the ground state 
as a function of rotation. 

The Hamiltonian for interacting bosons or fermions in a rotating optical lattice 
takes the second quantized form. 



where D is the dimension of the system, Viat (x) is taken to be sinusoidal. A, /x, and a are 
spin indices, and \E'ct(x), ^Pj^ (x) are field operators satisfying the (anti-) commutation 




, (x) (xO , 



(15) 
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X(x),^Mx')l = 5.^.5 (x-x'), 



(x),^., (x') 



0. 



(16) 



Assuming spin- independent contact interactions, a Hubbard Hamiltonian can be derived 
in a lowest-band, tight-binding approximation by expanding the field operators in the 
lowest-band Wannier functions, w(x — Xj) [?TJ. A basis better suited in the presence of 
rotation is given by the modified Wannier functions, 



Wj (x) = exp - 



^-"^^ J A (x') ■ dx'^ w {-x- Xj) , 



(17) 



where A (x') = f2 x x' [42j. This modification captures the effect of rotation at low 
rotation speeds |42j ; however, this approximation needs to be further modified at higher 
rotation speeds where the density of the Wannier functions is modified along with the 
phase. 

The result for spinless bosons is 



{id) ^ ^ 



(18) 



A similar derivation for spin- 1/2 fermions yields 



For both cases, i,j are site indices with {i,j) indicating a sum over only nearest 
neighbours, and is the radial position of the i'th site. The phase (pij is given by 
the expression, 



bij = — A (x ) ■ dx = -^[xiVj - Xjyi). 



(20) 



The parameters e, t, and U are the on-site zero-point energy, the hopping parameter 
and the on-site interaction energy, respectively, identical to those of the standard Bose- 
Hubbard Hamiltonian [H]. They are given by the expressions. 



d^x w*(x — Xj) 
d^x w*{x — Xj) 
U = g J d^x |w(x — Xj~ 



2M 
'2M 



V' + Vlat (X) J W(X-X,), 
V^ + Vlat (X) 1 W(X-X, 



(21) 
(22) 
(23) 
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where g is a. two-particle interaction strength. The parameters Ai and A2 arise due to 



the phase factor in equation 17 and are given by 



Ai = J dx w* (x — Xi) {x — Xif w (x — Xj) , (24) 

A2 = 2 J dx w*{x — Xi) {x — Xif w {x — Xi) , (25) 

where w {x — Xi) are one- dimensional Wannier functions. All of these parameters can 
be numerically evaluated for a lattice of specific shape, period, and depth [26]. If a 
harmonic trap is included, fi^ — * — where VLt is the trap frequency. 

In the rest of the paper, a lattice depth of Vq = IO-Er is assumed, and all parameters 
in the Hubbard Hamiltonians are numerically computed for this depth. 

3.1. Strongly-interacting bosons 

In this section, we describe hard-core bosons in a rotating optical lattice. The quasi- 
angular momentum of an eigenstate for a one-dimensional system is then defined and 
monitored as a function of rotation speed. The analytic results herein derived are 
consistent with a previous numerical treatment of a two-dimensional system [12]. 

In the limit of very strong interactions, t/ 00, a gas of bosons enters the hard-core 
boson regime. In an optical lattice, this regime can be characterized by using a number 
basis where the occupation number of each site is either / or / + 1, with / an integer. 
We can encode this fact in the Hamiltonian by formally changing the properties of the 
creation and annihilation operators. Here we consider the case of filling factors less than 
one, as the results for higher filling factors are qualitatively identical (see the appendix 
for the general case). This can be effected by stipulating on-site anti-commutation 
relations for the operators - i.e.. 



[ai,aj^i\_= ai,a'j_^^ =0. (26) 

In this case, the interaction term vanishes, although strong-interactions still implicitly 
exist in the Hamiltonian. For practical purposes, U = lOOt, where t is the tunneling 
energy, is large enough to enter this regime [42j . 
The Hamiltonian for a ring-lattice is given by 

H = -f2(t+ ^^1) + ^-c- 

i=i ^ ^ 

+ E - ^MQ' {R' + ^2)) a]a„ (27) 

where R is the radius of the ring. The parameter $ is given by 
M 

^ = ^{xiVj - XjVi), (28) 

which is a constant on a ring. This Hamiltonian can be diagonalized by first performing 
the Jordan- Wigner transformation [l3l HI] , 

6. = aj.e*"n=M««^ (29) 
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Figure 6. (Color online.) (a) The operator d[ adds one particle to the system, 
populating every site with equal probability but a phase that differs by 27r/8 between 
lattice sites. The state that is illustrated is jV') = ^1 |0); arrows indicate phase, 
(b) When this state is rotated by one lattice site, the resulting state is i?27r/8 IV')- 
Subtracting the phases at a site before and after rotation, we get = — 27r/8. The 
picture then illustrates the fact that \tp) is an eigenstate of i?27r/8 with eigenvalue 

g-i27r/8 



which transforms the mixed a operators into fully fermionic ones. The transformed 
Hamiltonian takes different forms for even and odd numbers of particles. Performing 
the canonical transformation, 

1 ^ 

^- = ^E^"^^^™'^%' (30) 



where fm = 2m for n odd, fm = 2m — 1 for n even, and n is the number of particles, 
results in 

N 

H=Y,Emdldm, (31) 

m=l 

= -2(^t+ ^MQ^Ai^ cos - Q<^^ + e - ^MQ^ (i?^ + A2) (32) 

The eigenstates are of the form dl^^ . . . d\^^ |0). 

A complete picture of this diagonalization requires an interpretation of the 
eigenstates in terms of quasi-angular momentum. The discrete rotation operator is 
defined via 

Rh]^^R-^ = h]^,, Rb^^R-' = i-irbl R\0) = \0). (33) 
This operator commutes with the Hamiltonian, and its action on an eigenstate, {ip) = 

R 1^) = e-^"(^-i+-+^-")/^ \^) . (34) 

Therefore, the quasi-angular momentum of this state is ^Iihh^L—^Ihul^ mod A^. A pictorial 
representation of this discussion is shown in figure 6 for a specific case with one particle. 



We are now in a position to use this formalism to investigate the ground state of 
the system. The ground state is determined by populating the lowest energy levels in 
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Figure 7. (Color online.) Ground state of three strongly-interacting bosons 
in a rotating 8-site ring lattice. Solid curves are guides for the eye. (a) At 
f2 = 0, the three lowest single-particle energies carry quasi-angular momentum 
m = — 1,0, 1. The three-particle ground state therefore has quasi-angular momentum 
m — (—1 + 0+1) mod 8 = 0. (c) At = 0.14ER/h, the three lowest single-particle 
energies carry quasi-angular momentum m = 0, 1,2. The three-particle ground state 
therefore has quasi-angular momentum m = (0 + 1 + 2) mod 8 = 3. 



the standard way for fermions. Equation 32 is the single-particle energy spectrum Em- 
As Q is increased, the minimum value of the cosine moves to the right, as illustrated 
in figure 7(a). The ground state is then determined by populating the lowest energy 
dln^rn eigeustates. In figures 7(b) and 7(c), this process is illustrated for the case of 
three particles in an eight-site ring. 

In general, if the minimum is near m = mo, we populate states mo, mo±l, . . . , mo± 
p for n = 2p + 1 particles in the system. The quasi-angular momentum of the state is 
then 

p 

m = (mo + /) mod = mo i2p + 1) mod = mo^ mod A^. (35) 
i=-p 

A similar argument holds for even numbers of particles. Thus, as rotation is increased 
from Q = 0, the quasi-angular momentum cycles through the values, 

m = jn mod A^ : j = 0, 1, 2, . . . . (36) 

One final note for the case of bosons. If there is an infinitesimal potential, 
V, that breaks the symmetry of the original Hamiltonian, Hq, down to a four- fold 
rotational symmetry, the behavior described in section 2.2 obtains. We treat the problem 
perturbatively. Tiny bandgaps in the single-particle energy spectrum, equation [32} open 
up at the new Brillouin zone boundaries. The eigenstates of Hq are independent of Q. 
Therefore, to first order in V, the energy level shifts are independent of Q, and the 
eigenstates remain unchanged. The effect of rotation is then exactly the same as in 
the non-perturbed case, and the quasi-angular momenta have the same behavior as a 
function of rotation. However, since the old quasi-angular momentum is no longer a 
good quantum number, we have to take 

m — i> m mod 4. (37) 

These results are consistent with the treatment in a previous paper |12] . 



Quasi-angular momentum in rotating optical lattices 



13 



3.2. N on- interacting fermions in a ring lattice 

Previous results focused on strongly-interacting bosons in rotating optical lattices 
[30l H2] . Here, we describe non-interacting fermions in ring lattices and compare to 
the boson case. The quasi-angular momentum of the ground state changes as rotation 
is increased. The allowed values of this quantity for fermions differ from those of the 
boson case. 

For non-interacting fermions, we drop the spin index, and the Hamiltonian for a 
ring lattice is identical to equation |27] except that the operators are fermionic. This 
Hamiltonian can be analytically diagonalized via the canonical transformation 

1 ^ 

^™ = 47E^""™''^%- (38) 



N . ^ 



The Hamiltonian is identical to that of equations 31 and 32 except that fm = 2m for 
both even and odd numbers of particles. The eigenstates, dl^_^ ■ ■ ■ d\^^ |0), carry quasi- 
angular momentum m = (mi + ■ ■ ■ + m„) modA^, where n is the number of particles. A 
discussion similar to the one given for bosons yields m as a function of rotation: 

m = nj mod A^; j = 0, 1, 2, . . . ; n odd, (39) 
m = n ( j -\ — ] mod A^; j = 0, 1, 2, . . . ; n even. (40) 



2^ 

As expected, the single-particle cases for both bosons and fermions are identical. 
It turns out that the cases of odd numbers of particles also coincide. However, there 
is an interesting distinction between the two cases when the number of particles in the 



system is even (compare equations 36 and 40). For instance, a system of two fermions 



in a four-site lattice cycles between quasi-angular momenta 1 and 3 whereas a system 
of two bosons cycles between values 2 and 4. This is one way in which non-interacting 
fermions and strongly-interacting bosons in one dimension differ despite the Jordan- 
Wigner transformation that maps between the two cases. 



3.3. Non-interacting fermions in a two-dimensional square lattice 

Differences between fermions and bosons remain when considering two-dimensional 
systems. Here we describe non-interacting fermions in a two-dimensional 4x4 square 
lattice. 

The Hamiltonian can be written as the sum of three terms, H = if 12 + -\- V, 
where H12 is the Hamiltonian of the outer twelve-site ring, is the Hamiltonian of the 
inner four-site ring, and V is the hopping between the rings, 

V = -J2(t + ^A,) ap)t«f e-*^ + h.c., (41) 

where the superscripts 4 and 12 indicate whether site j is on the outer 12-site or inner 
4-site ring. The Hamiltonians H4 and H12 can be separately diagonalized (apart from 
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Figure 8. (Color online.) Single-particle energy spectrum in a 4 x 4 square lattice for 
(a) zero rotation and zero hopping between inner and outer rings, (b) zero rotation and 
non-zero hopping, (c) = 0.05Efi/h and zero hopping, and (d) ft = 0.05Efi/h and 
non-zero hopping between inner and outer rings. The dashed and solid line are aids to 
the eye that indicate the band. In (a) and (c), the red dashed line is the spectrum for 
the inner 4-site lattice. 



the term in Hu involving rf) via the transformations, 

1 ^ 

= ^Y1 e"^2-'"^'/^af ) . (42) 
V ^ i=i 

In figure 8(a) is plotted the single-particle energy spectrum for the Hamiltonian in 
the absence of hopping between the rings {V = 0). When the hopping is "turned on" 
(figure 8(b)), level crossings become avoided crossing due to the mixing of states with 
the same four-fold quasi-angular momentum. In figure 8(c), the case for zero hopping 
between rings is plotted for Q = O.ObEn/h, whereas in figure 8(d), the hopping has 
turned on, and the avoided crossing behavior is observed. 

In figure 9 are plotted the quasi-angular momenta of the ground state as a function 
of rotation speed for 1, 2, 3, 4, and 5 non-interacting fermions in a sixteen-site lattice. 
Similar to the ring lattice case, the results for odd numbers of particles are identical to 
those for hard-core bosons, but differ for even numbers of particles. However, there is 
an extra distinction that does not occur in the ring lattice. The pattern of m values 
taken on for five particles (see figure 9(e)) differs markedly from that of one particle (see 
figure 9(a)), whereas for bosons these patterns are essentially identical [12] . 

4. Signatures of quasi-angular momentum in the momentum distribution 

While the quasi-angular momentum m is a good quantum number and a useful tool 
in investigating the symmetry properties of the ground state, it is not a quantity that 
can be directly measured in experiments. Instead, experiments are routinely performed 
in which time-of-fiight, far- field, density images of the gas contain directly information 
about the momentum distribution before expansion. As a first step towards detecting 
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Figure 9. (Color online.) Quasi-angular momentum as a function of rotation for (a) 
1, (b) 2, (c) 3, (d) 4, and (e) 5 non-interacting fcrmions in a sixteen site square lattice. 
These patterns repeat for larger rotation rates. 



quasi-angular momentum states, we look for a signature of quasi-angular momentum in 
the momentum distribution of a state. 

The momentum distribution of the many-body state, \ip), is given by 

n(k) = (^/'|^t(k)^(k)|V'), (43) 

where ^(k) is the Fourier transform of the field operator ^(x). Recalling that we expand 
the field operators as 

^(x) = ajWj{^) = J2 aje-'^'^w (x - x^) , (44) 

j j 

it can be shown that the momentum distribution takes the form 

n(k) = 5^^i;*(k;)«;(k,)e^''-(^'-^") a\a, |^) , (45) 
w 

where 

M 

k, = k + fi— x/xz. (46) 

This distribution is written in the momentum coordinates, k, corresponding to the 
rotating frame coordinates, x. Since the number density rotates in the lab frame 
coordinates, x^^, the momentum distribution also rotates in the lab frame momentum 
coordinates, k^. However, far-field pictures of the gas will be snapshots of the 
momentum distribution at the moment the trap and lattice are turned off, provided 
the switch-off time is fast enough. Thus, the momentum distributions presented here 
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are accurate representations of what will be measured in time-of-flight measurements, 
although they may be rotated relative to each other. 

Since quasi-angular momentum is a reflection of the symmetry of the system, 
considering single-particle states should be sufficient to capture the basics of how this 
quantum number affects the momentum distribution. We therefore concentrate on 
single-particle states only. In this case, n(k) is merely the square of the Fourier transform 
of the wavefunction. 

In order to build up an understanding of the momentum distribution, first drop 
the Wannier functions. In this case, n(k) is the Fourier transform of a sum of weighted 
delta functions. On a four site-lattice, the wavefunction is of the form 



(47) 



i=i 

where are the four corners of a square centered at x = 0. In figure 10 are plotted the 
momentum distributions for m = 0, 1, and 2. The m = is peaked at k = 27Tn/d, with 
n an integer and d the lattice spacing, whereas the other states vanish at these points. 
The spacing between peaks for both m = and m = 2 is Ak = 2TT/d, whereas for 
m = 1, 3 this spacing is Ak = 27r/-\/2c?; these can be calculated analytically. Thus, since 
the lattice spacing is known, there are clear measurable distinctions between different 
generic quasi-angular momentum states. Note that m = 1 and m = 3 are identical as 
they represent similar states with counter-propagating current patterns. 
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Figure 10. (Color online.) Fourier transform of generic (a) m — 0, (b) ni — 1,3, 
and (c) m = 2 quasi-angular momentum states on four sites. These momentum 
distributions are plotted over a range Air/d < kj.,ky < Aw/d; the first Brillouin zone is 
given by 7r/d < kx,ky < ir/d. (d) A schematic of the peak structure of the momentum 
distribution of four delta functions of the m = (blue) , m = 1,3 (green) , and m = 2 
(red) quasi-angular momentum states. The peak spacings for the m = and m — 2 
states are both 2t: /d, where d is the spacing between the delta functions, whereas the 
peak spacing for m = 1, 3 is 2-k /^/2d. 



Re-including the Wannier functions contributes to the momentum distribution an 
overall envelope of approximate width ^JvJE^/A, where V is the depth of the lattice. 
The modification k ^ k;, equation 46 moves the peak of this envelope away from k = 
as rotation is increased, though the center of the envelope is always located at k = 0. 
This is consistent with higher momenta being accessed for higher rotation rates. 

At this point, it is important to note that locating k = in the distribution is 
necessary for distinguishing the m = and m = 2 states. This can be done easily due 
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to the overall envelope in momentum space that is centered at k = 0. 

There are two additional considerations that come into play for larger lattices, both 
system-size effects. An overall envelope in position space determines the width of the 
peaks in momentum space. This envelope determines the size of the system, and if the 
width is roughly L, then the width of the peaks in momentum space scales as 1/L. If 
the lattice spacing is increased by a factor of n, then the peak-spacing in momentum 
space is decreased be a factor of 1/n. The general peak structure is not modified by 
these considerations. 

Since the peak structure is not modified by lattice-size effects, we can conclude that 
the distinctions between quasi-angular momentum states are sustained for more physical 
systems, which we turn to next. The ground state momentum distribution for a one- 
hundred-site lattice in the presence of a harmonic trap of frequency = O.lbEn/h is 
computed via imaginary time propagation. The trap ensures that the majority of the 
number density resides on the inner 36 sites. 

The results are summarized in figure 11, in which are plotted the momentum 
distributions for m = 0, 1, and, 2 states at rotation speeds of 0, O.lEu/h, and O.lAbEuh, 
increasing to the right. We can again make the observation that the m — state is the 
only one that is non-zero at exactly the reciprocal lattice vectors, i.e. k = ^k-r-|-^^ky. 
The peak-spacing is smaller for the m — 1 state than for the m = 2 state, distinguishing 
them from each other. However, the exact peak-spacing is unclear due to complex 
interference effects between sites at different radii. 

These differences arc washed out at higher rotation speeds. The m = 1, m = 2, and 
m = 3 (not pictured) states have very similar structures for large rotation. The peaks 
overlap so that it is very difficult to resolve them; the momentum distributions then 
all appear to be rings centered at reciprocal lattice vectors. They are distinguishable 
via the radius of the ring, as the radius is smallest for m = 1 and largest for m — 3, 
but as this is in order of increasing energy, the larger radius may just be an artifact of 
the larger width in the number density envelope at higher energies. Indeed, figures 11a, 
lib, and 11c are consistent with this conclusion, as the momentum peaks are narrowing 
due to the number density spreading out in real space. These results are consistent with 
those calculated via direct diagonalization of the Hubbard models previously discussed. 

5. Conclusion 

By analogy with quasi-momentum in translationally invariant periodic systems, the 
notion of quasi-angular momentum was introduced to label the eigenstates of a 
Hamiltonian that has a discrete rotational symmetry. It was shown that quasi-angular 
momentum is useful in analyzing the ground-state properties of quantum gases of bosons 
or fermions in rotating optical lattices. In particular, monitoring the quasi-angular 
momentum of the ground state as a function of rotation allowed us to identify transitions 
between different circulation values. 

We also presented a possible avenue by which the quasi-angular momentum of a 
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Figure 11. (Color online.) Momentum distributions computed via imaginary time 
propagation for one particle in a 100-site lattice in the presence of a trap of frequency 
fir = QlhEn/h. Plotted are m = ((a), (b), and (c)), m = 1 ((d), (e), and (f)), 
and m = 2 ((g), (h), and (i)) states for rotation speeds of = ((a), (d), and (g)), 
Q. = Q.lER/h ((b), (e), and (h)), and = QAA^Er/H ((c), (f), and (i)). At zero 
rotation, the distinctions between different quasi-angular momenta are sustained. As 
rotation is increased, it becomes harder to distinguish m — \ and m = 2 states. Due to 
the centrifugal term, higher rotation increases the width of the envelope of the number 
density, decreasing the width of the peaks in momentum space. Not shown is m = 3 
which has a structure similar to both m — \ and m = 2 for large rotation but has a 
ring-radius larger than both of them. 



state can be experimentally determined. We identified characteristics in the momentum 
distribution distinguishing between different quasi-angular momentum states at low 
rotation speeds, such as the existence of a peak at reciprocal lattice vectors for an m = 
state only, the peak-spacing, and the overall structure of the momentum distribution. 

There are still open questions as to how the momentum distributions will change 
when the lattice size or number of particles is increased. The effects of statistics will 
be of fundamental importance for many-particle systems; there is a question of how 
well the peaks can be resolved for larger lattices, since peak-spacing decreases with 
increasing lattice size; and interference effects between adjacent sites and between rings 
is not rigorously treatable. However, there will always be signatures of the quasi-angular 
momentum in the momentum distribution, as it is connected with the phase information 
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of the ground state, which in turn influences the momenta of the system. In addition, 
a zero quasi-angular momentum state can be distinguished from a non-zero one, which 
allows the experimenter to verify that vorticity has entered the system. 
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Appendix 

The case of fiUings between I and I + 1 can be encoded in the Hamiltonian by formally 
changing the properties of the creation and annihilation operators. As mentioned in 
the body of the paper, for fillings less than one, the creation and annihilation operators 
satisfy the following relations: 

j' 



[ai,aj] = ai,a] = 0. (A.l) 

If the filling fraction is between / and / + 1, the following changes are made. First, the 
site number operator ajaj is replaced by UjUj + /, reflecting the fact that each site is 
filled with between I and I + 1 particles. Secondly, recalling that the operators were 
originally bosonic, Oj \l) — -s/T+T |i -l- 1) and aj\l + 1) — y/l + 1 \l), one can see that 
the hopping parameter t is scaled by a factor of Z -|- 1. Finally, the interaction term gets 
modified according to the replacement Ojttj OjUj + I: 

N N 

a]aj {a]aj - l) ^ 2/ ^ a]aj + NI{1- 1) (A.2) 
j=i i=i 
The Hamiltonian in two-state approximation is then given by 

^ 1(1-1) ^ 

H ^ {Ul - - A*^^ + UN ^ ^ ' - (Z + 1) ^ te"l'a]^^aj + h.c. (A.3) 

3=1 j=i 
where the chemical potential /i has been introduced so that we can work in the grand- 
canonical ensemble. 

Upon applying the Jordan- Wigner transformation, 

6j = a^•e''^S'fc=l«^«^ (A.4) 

the Hamiltonian becomes 

N , /, , N N-1 



1(1-1) 

H ^ {Ul - b]bj - f^Nl + UN ^ ^ ' -{l + l)Y^ te^l+ibj + h.c. 

j=l 3=1 

- (-1)"+' (/ + 1) te%{hN + h.c. (A.5) 
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where n is the number of particles. Now that the Hamiltonian is fully fermionic, we can 
perform a canonical transformation to diagonalize it, 

1 ^ 

rfm = ^5^e^"^-'^/%, (A.6) 



where fm = 2m for n odd and fm = 2m — 1 for n even. The resulting Hamiltonian is 

N 



H = -fiNl + (/ - 1) + ^ Emdld„ 



m=l 



E^ = -^-2{l + l)tcos['^-(l)] +Ul. (A.7) 



Setting / and /x to zero, we realize the Hamiltonian derived in Sec. 3. 

As an aside, this calculation actually yields more than just the description of the 
quasi-angular momentum of a rotating ring lattice. In the limit of large A^, this system 
is indistinguishable from a linear lattice. Taking t real, the Mott-insulator/superfluid 
phase diagram in two-state approximation can be derived. Allowing the number of 
particles to vary, it is apparent that the number of particles is determined by the number 
of m's for which Em < 0. Varying the chemical potential /x, every time there exists an m 
such that Em = 0, there is a phase boundary between ground states of different numbers 
of particles. Thus, the phase boundaries are given by the expression, 

^ - -2(/ + l)^cosf^V^- (A-8) 



U ' ' U \ N ^ 

In particular, it is possible to identify regions of integer filling and non-integer filling, 
the boundaries between which are given by 

^ = ±2(1 + 1)^ + 1 (A.9) 

for a large number of sites. These boundaries yield an approximation to the Mott- 
insulator superfiuid phase diagram (figure Al). They have been derived in a slightly 
modified form in using a semiclassical approach. We note that our method can 
be easily extended to the case of superlattices, generating for instance Mott phases of 
half-filling in the phase diagram. 
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